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ABSTRACT 

We study, for the first time, how shear and angular momentum modify typical parameters of 
the spherical collapse model, in dark energy dominated universes. In particular, we study the 
linear density threshold for collapse <5 C and the virial overdensity Ay, for several dark-energy 
models and its influence on t he cumulative mas s function. The equations of the spherical 
collapse are those obtained in iPace et all (1201 Oh . who used the fully nonlinear differential 
equation for the evolution of the density contrast derived from Newtonian hydrodynamics, 
and assumed that dark energy is present only at the background level. With the introduction 
of the shear and rotation terms, the parameters of the spherical collapse model are now mass- 
dependant. The results of the paper show, as expected, that the new terms considered in the 
spherical collapse model oppose the collapse of perturbations on galactic scale giving rise to 
higher values of the linear overdensity parameter with resp ect to the non-rotating case. We find 
a similar effect also for the virial overdensity parameter. For what concerns the mass function, 
we find that its high mass tail is suppressed, while the low mass tail is slightly affected except 
in some cases, e.g. the Chaplygin gas case. 
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1 INTRODUCTION 

Till a decade ago, Universe was considered composed mainly 
by dark matter (DM) and characterized by a decelerating expan- 
, sion. An important and surprising result coming from observa- 
tional cosmology is the fact that high redshif t supernovae are 
less bright than expecta tions dRiess et al] 1 19981 : IPerlmutter et al.l 

1 19991 : iTonrv et al.l l2003h . This finding has been interpreted as 
an acceleration in the expansion of the universe and that this 
acceleration is recent (Riessetal. [l998: Perl mutter et al.l [l999l : 

iKnopetalJ [20031 : iRiess et alj 120041 : 1 Astier et alj 120061) . This re- 
sult has been confirm ed by independent observations: the baryon 
acoustic oscillation dTegmark et al.l [20041 : lEisenstein et alj [20051 : 

IPercival et al.l2O10l) . the galaxy-galaxy correlation function, giving 
important infos on the spatial distribu tion of large-scale structure 
(Tegm ark et al.l20 04: Cole et al. 2005), the angular spectrum of the 
CMBR temperature fluctuations (Komatsu et alj201l | : |Larson et al. 

l201lh . the Integ rated Sachs-Wolfe (ISW) effect dHo et . al.ll2008f> 7 
globular clusters dKrauss & ChaboveifeOoHlDotter et al.l201 lh . old 
high redshift galaxies (A l caniz et al. 2003, 2005), and galaxy clus- 
ters jHaiman et al.ll200l1 : lAllen et al.ll2004 120081: iBasilakos et alj 

|2010|) and weak lensing jHoekstra et alj|2006l : Ijarvis et al.l l2006). 
The quoted accelerated expansion cannot be obtained in uni- 
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verses containing just matter homogeneously and isotropically dis- 
tributed, while it can be obtained if the low-z universe is filled 
in with a fluid with negative pressure, the so called dark energy 
(DE), with equation-of state-parameter w < —1/3. It is possi- 
ble to have accelerated universes without dark energy if one dis- 
cards the homogeneity hypothesis on large scales (e.g., Lemaitre- 
Tolman-Bondi ( LTB) universes), using the back-reaction approach 
to dark energy dKolb et alj |2006). gravitationally induced particle 
creation (Lima et al. 2009, 2 010) or even mo dified models of grav- 
ity, like the f(R) models dAmendola et alj 120071). f(T) model s 
(Bengochea & Ferrarol2009t) . or the brane models jDeffavetl200ll) . 

The nature of dark energy is not understood to date, and this 
explains why in the past decade a large number of models for the 
origin and time evolution of DE have been proposed. In the ACDM 
model the DE is connected to the energy of vacuum (cosmological 
constant) and the equation of state of DE, in this case, is simply 
w — —1. An extension of this model is obtained considering a 
scalar field with no or weak interaction with the matter component 
(quintessence models), and also phantom models, K-essence, or al- 
ternatively Chaplygin gas and Casimir effect. 

Cosmologists generally believe that structures in the universe 
grew via gravitational instability through the growth and collapse of 
primeval density perturbations originated in an inflationary phase 
of early Universe from quantum, Gaussian distributed, fluctuations 
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I Guth & Pill982l ; lHawkinelll982l : IStarobinskvlll982l ; lBardeen et al.l 



Il986h 



The presence of DE changes the rate of formation and growth 
of collapsed structures and large-scale structure, and consequently 
the distribution in size, in time and space of galaxies, quasars, su- 
pernovae, since they reside in collapsed structures. Moreover, DE, 
increasing the expansion rate, slows down the collapse of overdense 
structure and its space-time fluctuations (if DE is not the cosmolog- 
ical constant) will give rise to DE haloes dCreminelli et al.ll2010l) 
which will influence dark halos formation. 

In principle, all the statistical information concerning the dis- 
tribution of the dark matter and DE is contained in the probability 
distribution functions (PDFs) for the velocity, v, and mass density 
fluctuation fields, 8. The determination of the final moments of 
the quoted distribution from the starting ones require to know the 
exact dynamics ruling the evolution of the underlying field. Unfor- 
tunately the exact solution to the dynamical equations is known 
only in the linear regime of evolution of the quoted fields, and only 
approximate solutions are known for non- linear stages. A popular 
analytical approach to study the non-linear evolution of perturba- 
tions of DM and DE is the sp herical collapse model (SCM) intro- 
duced in the seminal paper of iGunn & Gotq dl972» extended and 
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1986, 1989'; IZaroubi & Hoffman! 
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issue of the role of velocity shear in the gravitational collapse. 

The model describes how a spherical symmetric overdensit)0 
decouples from the Hubble flow, slows down, turns around and col- 
lapse. In the model the overdensity is divided into bound mass, each 
one expanding with the Hubble flow from an initial comoving ra- 
dius Xi to a maximum one x m (usually named turn-around radius, 
x ta ), and then collapse. Non-linear processes convert the kinetic 
energy of collapse into random motions, giving rise to a "virial- 
ized" structure. 

The SCM proposed by IGunn & GotH d 19721) does not 
contain non-radial motions and angular momentum. The way 
to introduce angular momentum in the SCM, and its con- 
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purely radial self-similar collapse by including non -radial motions 
arising from secondary perturbations in the halo. [Williams et al.l 



(2004) used the same model to show how angular momentum flat- 
tens the inner profile of haloes. The effect of non-radial motions 
on the mass profile in a SCM was studied by iGurevich & Zvbinl 
dl988allbh 



White & Za ritskv! l l 1992b applied a torque to particles in the 



shells during the initial expansion phase, and in order to preserve 
the spherical symmetry, assumed that the different particles ac- 
quire the same angul ar momentum but in independent randomly 
oriented directions. In Nusserl d200ll) . particles acquire an angular 
momentum at turn around, while before they move on radial or- 



1 A slightly overdense sphere, embedded in the universe, is a useful non- 
linear model, as it behaves exactly as a closed sub-universe because of 
Birkhoff 's theorem. 



bits. Again, in order spherical symmetry is preserved and the angu- 
lar momentum of each particle is conserved, particles have angular 
momenta distributed in random directions such that the mean angu- 
lar momentum at any point in sp ace is zero. More over, angular mo- 
mentum is oc \/GM(< r«)rjf| per unit mass IWhite & Zaritskvl 
ll992l ; ISikivie et al.lll997l). so no additional physical scale is intro- 
duced. Other studiesTHiotelisl2002l ; lLe Delliou & Henrikseij2003l: 
lAscasibar et alj|2004h introduced angular momentum in the SCM 
in a similar way of the previous cited authors, and studied the struc- 
ture of DM density pr ofiles, reaching similar conclusions to that of 
IWilliamsetai1d2004l) . 

The SCM in the framework of d a rk en e rgy cosmolo- 
gies was studied by iMota & van de Bruckl d2004l) ; lAbramo et al.l 
d2007l); iBasilakos et alj d2010r i: IPace et alj \201(E In particular in 
iPace etalJd2010h were derived the equations for the SCM under the 
assumption that only DM c an form clumps and that DE is present as 
a background fluid (see alsolFosalba & Gazta fiaga 1998; Oht a et ail 
l200il2004l ; lMota & van de Bruckll2004l;lAbramo et al.ll2007l) . The 
evolution of the overdensity S is given by (Pace et al. 2010, 2012): 



8 + 2H5 



4 5 2 



31 + 5 



-4nGp8(l + 5)- 
(1 + 5){a 2 - ui 2 ) = 



(1) 



aija 1 - 1 and the rotation term uj — ooijio 1 - 



the shear term a 2 
are connected to the shear tensor, which is a symmetric traceless 
tensor, while the rotation is antisymmetric. They are given by: 



— — 



dx 1 
du 3 
dx l 



~dxi 



(2) 



(3) 



where 8 = Vs • u is the expansion. 

Recalling that 5 = p/~p — 1 = (a/i?) 3 — 1 (a is the scale factor 
and R the radius of the perturbation), and inserting it into Eq.Q] it 
is eas y to check that the evolution equation for 5 reduces to the 
SCM dFosalba & GaztaiiagJl998l ; lEngineer et al.l2000l : IOhta et al.1 
2003) 



d 2 R 



— = AnGpR - l/3(<7 



- U )R = -—. 



-l/3(a 2 



-u )R , 
(4) 

comparable with the usual expression for the SCM with angular 
mom entum (e.g. |Peeblesll993l ; lNusser 2001 ; Zukin & Bertschinger 
l2010h : 



d 2 R 
dt 2 



GM 



+ 



R 2 M 2 i? 3 



GM 4- 4 n 2 R 



(5) 



where in the last expression we have used the momentum of inertia 
of a sphere, / = 2/5MR 2 . 

The previous argument shows tha t vorticity, u>, is strictly con- 
nected to angular velocity, Q (see also lCherninl ([l993), for a com- 
plete treatment of the interrelation of vorticity and angular momen- 
tum in galaxies). 

One assumption generally used when solving the SCM equa- 
tions for the density contrast 8 (Eq. Q3 is to neglect the shear, a, 
and the rotation oj. While the first assumption is correct, since for 
a sphere the shear tensor vanishes, the rotation term, or angular 
momentum is not negligible. In fact, if we consider the ratio of the 
rotational term and the gravitational one in Eq.|5]we get M 3 RG that 



r* is the maximum radius of oscillation of a particle 
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for a spiral galaxy like the Milky Way, with L a ; 2.5 x 10 74 gem 2 /s 
dRvden & Gunrj(l987l ; ICatelan & Theunsll 19961) , and radius 15 kpc 
is of the order of 0.4, showing, as well known, that the rotation is 
not negligible in the case of galaxy sized perturbations. The quoted 
ratio is larger for smaller size perturbations (dwarf galaxies size 
perturbations) and smaller for larger size perturbations (for clusters 
of galaxies the ratio is of the order of 10~ B ). The value of angular 
momentum, L, or si milarly Q., can be obtained and added to the 
SCM as described in lDelPopoldJ2009l) or as described previously, 
assi gning an angular momentum ex \/ GM(< r,)r* at turn-a round 
(e.g lwhite & Zaritskvll 19921 : Isikivie etZlll99llNusse^l200ll) . 

As previously stressed, the non-trivial role of angular momen- 
tum in the SCM has been pointed out in a noteworthy number of pa- 
pers studying structure formation in DM dominated universes (see 
also DelPopold 120091 : IZukin & Bertschingej|2010l : ICupani et all 
I20T 1|). In a previous letter. iDel Popolo et all d201 3) studied the ef- 
fect of the term a 2 — oj 2 on the SCM parameters (S c and Ay) for 
the Einstein-de Sitter (EdS) and ACDM models, but it has never 
taken into consideration in the SCM in DE cosmologies. 

In the present paper, we shall study how the typical parame- 
ters of the SCM (in Universes dominated by DE), namely the linear 
density threshold for collapse S c and the virial overdensity Ay, are 
changed by a non-zero a and uj terms. In fact, any e xtension of the 
SCM should take into accou nt the effects of shear (Engine er et al.l 
l2000l : lDel Popolo et alj2012h since shear induces contraction while 
vorticity induces expansion as expected from a centrifugal effect. 
We also study how angular momentum and shear influence the cu- 
mulative mass function. 

The paper is organized as follows. In Sect. [2] we summa- 
rize the model used to obtain 5 C and the virial overdensity Ay. 
In Sect. [3] we describe the results, and Sect.|4]is devoted to conclu- 
sions. 



2 SUMMARY OF THE MODEL 

The evolution equations of 8 in the non-linear regime has been 
obtained and used in the framework of the spherical and ellip - 
soidal collaps e, and st ructure formation by Bernardeau (1994); 
IPadmanabhai] dl996h : bhta et al.l d200j |2004 ): lAbramo et al.1 
j2007l) . As a first step we assume that the fluid satisfies the 
equation-of-state P = wpc 2 . In addition, w e also consider the 
Neo-newtonian expressions dLimaetal.lll997l) for the continuity, 
the Euler equations, and the relativistic Poisson equation, namely: 



dp „ , ^ P 



Vr • V 







dv 

v 2 $ 



+ (v ■ Vf)v + = 



4ttG p + 



3P 



= 



(6) 



(7) 



(8) 



where v is the velocity in three-space, $ is the Newtonian gravita- 
tional potential and r is the physical coordinate. 

The continuity equation for the mean background density can 
be written in the form 



P + 3H [p + 







(9) 



where p — 3 " fJ 8 ^g" id is the background mass density of all contri- 
butions to the cosmic fluid, and fWid is its density parameter. 

Using comoving coordinates x — f/a, the perturbations equa- 



tions can be written as 

8 + {l + w){l + S)V s -u = 

: 



du nrr ^ ,^ „ , _, 1 „ 
— + 2Hu + (u ■ V 3 )u + 
at cr 



V|0 - 4ttG(1 + 3w)a 2 p8 







(10) 
(11) 
(12) 



where H(a) is the Hubble function and u(x, i) is the comoving pe- 
culiar velocity. Combining the previous equations, we get the non- 
linear evolution equation 



5 + [2H - 



5- 



4 + 3w 5 2 



(13) 



1 + w J ' 3(1 + w) 1 + 8 

4ttGp(1 + w)(l + 3w)8(l + 6)- 

(l + w){l + 5)(a 2 -lj 2 ) = 0. 

which is a generalization of Eq. 7 of lAbramo et al.l d2007l) to the 
case of a non-spherical configuration of a rotating fluid. 
In the case of dust (w — 0), Eq.|13|reads 

4 S 2 

u 47rGp5(l + 5)- 

H y ' (14) 



8 + 2HS 



31+8 ' v ' 

(1 + 8)(o 2 - uj 2 ) = . 

which is Eq. 41 of lOhta et alj d2003l) . 

In terms of the scale factor, a, the nonlinear equation driven 
the evolution of the overdensity contrast can be rewritten as: 

„ [3 E'\ , 4 S' 2 3 n m , 

(1 + 8){a 2 - to 2 ) = , 



a?H 2 {a) 



(15) 



where Q m , is the density parameter of the DM at a = 1, E(a) is 
given by 



E(a) 



H f- +^Q,og{a) 



where g(a) is 



g(a) = exp ^~ 3 ^ 



a l + w(a') , , 
— — - da 



(16) 



(17) 



Note that in Eq.[T5]we corrected a typo present in Eq. 17 of 
IPace etalJd2oToh . 

In order to calculate the threshold for the collaps e # c and th e 
virial overdensity, Ay, of the SCM, we follow IPace et alj d2010l) . 
We look for an initial density contrast such that the 8 solving the 
non-linear equation diverges at the chosen collapse time. Once the 
initial overdensity is found, we use this value as an initial condition 
in the linearised equation 

\a E J 2 a 5 E 2 

to get 8 C . 

The initial conditions to solve the second-order differential 
equations are Si (got as previously described) and the initial rate 
of evolution, 8{ is calculated as follows. We assume that at early 
times the solution is a power law, therefore we can write Si = Ba n , 
where n is in general of the order unity. The velocity is defined as 
the derivative with respect to the scale factor of the initial overden- 
sity, hence we can write it as 8[ = n8i/a. 

The virial overdensity is obtained, as in IPace et alj d2010h . by 
using the definition Ay = log (6 n \ + 1) = (,(x/y) 3 , where x = 
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a fata, is the normalised scale factor and y is the radius of the sphere 
normalised to its value at the tum-around. 

The turn-around scale factor is obtained by solving Eq.ll5land 
determining the quantity log(<5 n i + l)/a 3 . The virial overdensity at 
tum-around £, is obtained by integrating Eq.[T5]up to a ta and add 
the result to unity. 

In order to integrate Eq. 15, we should explicit the a 2 — uj 2 
te rm. The calculation of th e term a 2 — uj 2 is explained in detail 
in iDel Popolo et alj d2012h . Here we simply summarize how we 
evaluated it. 

We first define with a the dimensionless ratio between the ro- 
tational and the gravitational term in Eq. ((5}: 



having different values for different scales, as already reported. 

We may calculate the same ratio between the gravitational and 
the extra term appearing in Eq. 15 thereby writing the term extend- 
ing the standard SCM as 

a 2 - uj 2 3 afi m ,o x 

-fff- = -2^^' (20) 

as recently discussed in the literature (Del Popolo et al 
In order to obtain a va lue for 8 C similar to the one obtained by 
IShefh & Tormenl dl999h . we set q = 0.05 for galactic masses 
(M ~ 10 11 Mq) corresponding to a rotational velocity of v r w 
250 km/ s and scaled it linearly towards higher masses and low 
velocities by assuming a rotational velocity of nearly 10 km/s for 
galaxy cluster size objects (M w 10 15 Mq). 

Given the above ansatz, the nonlinear equation that has to be 
satisfied by the overdensity contrast S for a large class of noninter- 
acting dark energy models reads: 

\a E J 31 + 5 2 a 5 E 2 (a) v ; 

(21) 

It should be remar ked that the DE models used are the same 
of IPace et al.1 J2010h . namely the ACDM model, the quintessence 
models, phantom models and topologi cal defects, and th e Chaply- 
gin gas and Casimir effect. We refer to IPace etal.1 ( bold) for more 
details on the models used. 



3 RESULTS 

In this section we present results for the two main quantities derived 
in the framework of the SCM, in particular the linear overdensity 
parameter 5 C and the virial overdensity Ay. We assume as refer- 
ence model the ACDM model, with the following cosmological 
parameters: Q m = 0.274, f2 dc = 0.726 and h = 0.7. 

First of all, in Fig.[T] we show the linear overdensity parameter 
5 C for the EdS model (upper panel) and ACDM model (lower panel) 
as a function of the redshift and of the mass of t he collapsing ob- 
ject. E ven if a projection of this plot was shown in Del Popolo et al.1 
we here propose the surface spanned by 8 C to briefly sum- 
marize how this quantity depends on the mass and on the redshift. 

An important result that it is worth to notice, is the time be- 
haviour of S c . We observe that the contribution of the term a 2 — uJ 2 
is maximum at z — and it decreases with increasing redshift 
till S c reaches an approximately constant value, generally higher 
or equal to the value for the standard SCM, according to the mass 
range considered. This is expected when we compare the nonlin- 
ear term with the gravitational term. The net result is that of giving 




Figure 1. Linear overdensity parameter <5 C as a function of mass and redshift 
for the EdS model (upper panel) and ACDM model (lower panel). 

as source term a model with an effective matter density (1 — a) 
times smaller than the real matter density Q m (only in the non- 
linear regime though). This can be interpreted as an additional term 
counteracting the collapse even at high redshifts, making therefore 
S c higher than the standard value. 

In Fig. [2] instead we present results for the different dark en- 
ergy models considered in this work. On the left panel we show 
results for the linear overdensity parameter S c while on the right 
panel we show the expected values for the virial overdensity Ay. 
Because of the consideration expressed above regarding the EdS 
model, we will consider as reference the standard ACDM model 
(with a = 0) and to maximize the effect of the non linear term, all 
the figures show results for galactic masses (M « 10 11 M Q /h). 

For the first class of models (INV1, INV2, 2EXP, AS, CNR, 
CPL, SUGRA) we see that the models are in general very similar to 
each other and they slightly differ from each other, with differences 
at most of the order of 4.5% for the AS model. The DSTV1, CNR 
and 2EXP models are basically indistinguishable from the ACDM 
model. We interpret this result as d ue to the equation of state of 
the models considered. As shown in IPace et al l <20ich . the INV1 
model has an equation of state basically constant over the whole 
cosmic history, but its present value is quite different from all the 
other models, being wo ~ —0.4, while for all the other models is 
-1 «S w < -0.8. 

Comparing our present r esults with the ones of the upper left panel 
in Fig. 4 of IPace et al. we see that the behaviour of the 
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Figure 2. The left panels show the time evolution of the linear overdensity S c (z), the right panels the time evolution for the virial overdensity Av(z) for 
the different classes of models. In all panels, the ACDM solution (black solid curve) is the reference model. All the curves assume a galactic mass for the 
collapsing sphere. The upper panels present the quintessence models: the INV1 (INV2) model is shown with the light-green dashed (dark-green short-dashed) 
curve, the 2EXP model with the blue dotted curve, the AS model with the cyan dot-dashed curve, the CPL (CNR) model with the red dot-dotted (orange 
dot-short-dashed) curve and finally the SUGRA model with the violet dot-dot-dashed curve. The middle panels show the Casimir effect (brown dotted curve) 
and the (generalized) Chaplygin gas with the (turquoise short-dashed) magenta dashed curve. Finally the lower panels report the solution for the models with 
constant equation-of-state parameter: the dark-green short-dashed curve stands for the non-flat ACDM model, the light-green dashed curve for the model with 
w = —2/3, the blue dotted curve represents the model with w = —4/3, the cyan dot-dashed curve the model with w = —1.5, the orange dot-short-dashed 
curve the model with w = —5/3, the red dot-dotted curve the model with w = —2 and finally the violet dot-dot-dashed curve shows the model with w = — 3. 
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models is very similar. The inclusion of the non-linear term just 
changes the values of the linear overdensity parameter, but not the 
respective ratios with the ACDM model. 

For the second group of models (Casimir and (g eneralized) 
Chaplygin gas) we obtain very different results from IPace et all 
(2010). While there only the generalised Chaplygin gas was sub- 
stantially different from the ACDM model, now all the models here 
considered differ much from the reference model. This shows how 
the non-linear additional term is very sensitive to the equation of 
state considered. Moreover, none of the models recovers the ex- 
tended ACDM model at high redshifts. 

The bottom panel in the left column is devoted to the phan- 
tom models (w < —1) and to a non-flat ACDM model. All the 
models present very similar results and small differences appear at 
small redshifts (z < 1), and for redshifts z > 2, all the models 
are identical. Models differing most at z ~ are the models with 
w — —2/3 and w = — 3 having values for 8 C respectively lower 
and hi gher than the AC DM model. This is in agreement with re- 
sults of iPaceetal]b010h . showing once again that a super-negative 
equation-of-state affects only slightly the structure formation pro- 
cess. In particular, the more negative it is, the higher is the linear 
overdensity parameter. We also notice that a small amount of cur- 
vature, does not influence our results significantly. 

Comparing these results, with the results from Fig. 4 in 
IPace et all feoiO), we can appreciate the interplay between the term 
a 2 — cu 2 and a dynamical dark energy equation of state. 
All the models studied with a time-varying dark energy equation- 
of-state parameter show that the collapse, even if retarded by the 
inclusion of the shear and rotation, is easier as compared with the 
ACDM model. In this case, with easier we mean that the values for 
the extended S c (z) are smaller than for the reference model. This 
is expected and it has the same explanation as for the usual case. 
Since at early times the amount of dark energy is higher, we need 
structures to grow faster in order to observe cosmic structures to- 
day. This is plausible, since the linear overdensity parameter repre- 
sents the time evolution of the initial overdensity, whose evolution 
is dictated by the growth rate that is described by the same differ- 
ential equation. In other words, since at early times the amount of 
dark energy is higher, we need lower values of S c to have objects 
collapsing. This is analogous to the case of the linear growth factor, 
since the equation to be solved is the same. 

Opposite is instead the behaviour for the phantom models, in which 
case we notice that with a more negative equation of state, the col- 
lapse is retarded more severely. This is in agreement with Fig. 4 of 
IPace et all d2010l ). where phantom models had an higher <5 C : since 
the expansion goes so fast, the collapse is strongly suppressed and 
with respect to the reference models, higher and higher initial over- 
densities are required in order to have collapsed objects today. 

On the right panels, we show results for the quantity Ay. Also 
in this case we limit ourselves to galactic masses. We notice that 
for the quintessence models in general the virial overdensity shows 
higher values than for the extended ACDM model, except for the 
INV2 model. This is opposite to what found for the standard case, 
where all the models had smaller values than the ACDM model. 
Also in this case none of the models approximates the extended 
ACDM model at high redshifts. This is not the result of the addi- 
tional non-linear term only, but also of the influence of the dark 
energy equation of state, consistently with results from lPace et al.l 

J20ior> . 

For the phantom models, results are very similar to the usual 
case. Virial overdensity parameter is higher than the extended 
ACDM one if w < —1 and lower for the model with w = —2/3, 



Table 1. Table with the power spectrum normalization for the different dark 
energy models (K is the curvature parameter). 



Model 




^8 


ACDM, K= 





0.776 


ACDM, 





0.793 


INV1 




0.428 


INV2 




0.707 


2EXP 




0.739 


AS 




0.319 


CNR 




0.732 


CPL 




0.444 


SUGRA 




0.578 


Chaplygin gas 


0.066 


g. Chaplygin j 


jas 


0.133 


Casimir 




0.420 


Phantom (w = - 


-2/3) 


0.674 


Phantom (w = - 


-4/3) 


0.834 


Phantom (w = - 


-1.5) 


0.854 


Phantom (w = - 


-5/3) 


0.870 


Phantom (w = 


-2) 


0.894 


Phantom (w = 


-3) 


0.936 



in agreement wit h the linear overdensity parameter. As shown in 
IPace et at higher redshifts, all the phantom models re- 

duce to the ACDM model. 

Differences in the linear overdensity parameter reflect in the 
differential mass fu nction. In this work w e decide to use th e 
parametric form by Isheth & Tormenl dl999h : Isheth et al.l fcOOlh : 
ISheth & Tormenl d2002l) . We consider three different redshifts, 
namely z = 0, 0.5, 1, 

At this point, it is w orth noticing that in the first paper 
of ISheth& Tormenl i ll 9991) the mass function was calculated as 
a fit to numerical simu lations. Later on, Isheth et aL I J200lh and 
ISheth & Tormenl d2002l) showed that the effects of non-sphericity 
(shear and tides) intro duce a mass depe ndence in the collapse 
threshold (see Eq. 4 in Isheth et al. I J200lh and following discus- 
sion). By using this threshold as the barrier in the excursion set 
approach one gets a mass function in good agreement with simula- 
tions (see Eq. 5 a nd the discussi on in the last part of Sec t. 2.2. of 
ISheth et alj feOOlh . and moreover lshefh & Tormenl d2002l) ). 

For the ACDM model we choose as power spectrum normal- 
ization the value as = 0.776. Since we want that perturbations at 
the CMB epoch are the same for all the models, we normalize the 
dark energy model according to the formula 

£ ) +,acdm(cicmb) .... 

0"8,DE = 0"S,ACDM — =: ; r — , (22) 

-L-*+,DE(aCMB J 

where D+(acMB) is the growth factor at the CMB epoch. 

In Table Q] we display the different normalizations used for 
the models considered in this work (K is the curvature parameter). 

In Fig. [3] we compare the differential mass function for the 
ACDM model in the standard and in the extended SCM. 
Analysing the three curves, we can appreciate the effect of the term 
a 2 — cu 2 on the mass function. On small masses, the mass func- 
tion is largely independent of the cosmological model, but it de- 
pends strongly on S c . Since in the extended SCM this is higher, we 
observe an increase in the number of objects at galactic scale at 
z — 0, to decrease at higher redshifts where the contribution of the 
non-linear term decreases. We observe a general decrement in the 
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Figure 3. Ratio between the differential mass function of the extended and 
standard ACDM model. The curves represent three different redshifts: z = 
(dotted red curve), z = 0.5 (short-dashed blue curve), z = 1 (dashed 
green curve). 



number of objects at high masses (up to M ~ 10 14 Mq/K) to in- 
crease again to unity for masses of the order of 10 lj Mq /h. This is 
explained by the fact that at such masses, the linear overdensity pa- 
rameter are practically the same, therefore the mass function must 
not change. 

In this concern, it should be recalled that shear and rotation 
have the maximum effect on 5 C , at galactic scale (see Fig.Q~|l. How- 
ever, in the calculation of the ratio between the differential mass 
functions (Fig.[3), beside S c we need to take into account the factor 
a(M), the r.m.s. of mass overdensity. Now, recalling the Sheth & 
Tormen multiplicity function 



fsT = A 



1 + 



o(Mf 
52 a 



<t{M) 



(23) 



with A = 0.3222, a = 0.707, andp = 0.3, we have that at galactic 
scale the dominant term in Eq.[23]is the term — , consequently the 

ratio /sT, C xtcndcd//sT,ACDM — 5 c ,extended/5c,ACDM • At larger 

masses, the effect of rotation and shear diminishes with the con- 
sequence that 5c, extended — S c , acdm and that the a(M) term is 
fixing the value of the /sT,extended//si\ACDM ratio. 

At this point it is necessary a deeper discussion of the results 
shown in Fig. [3] The ST mass function generalises the Press & 
Schechter (PS) (Press & Schechtenll974l) mass function to include 
the effects of shear and tidal forces with respect to the simpler 
spherical collapse model. In doing so it is necessary to consider the 
ellipsoidal collapse model and the corresponding linear overden- 
sity parameter 5 C c- Differently from the spherical collapse model, 
now Sec is not only a function of ti me, but of mass t oo and the 
relation bet ween 5 BC an d S c (5 SC in ISheth et alj feOOlh ) is given 
by Eq. 4 in IShethetal.1 bOOlh . The moving barrier for the ran- 
dom walk is now set equal to S ec (M, z) and a good approximation 
to it is then given by their Eq. 5. The ST mass function fits quite 
well the results of N-body simulations, but as shown in Fig. [3] for 
masses M « 1O 14 M0 /h our predictions, including rotation on top 
of the ellipsoidal collapse, predict approximately 40% less objects 
at 2 = 1 than the standard ACDM model. This would be easily 
checked with a big enough spanned volume. Moreover one could 
identify ellipticity and rotation if halos acquire angular momentum 
by misalignment with the surrounding tidal field. Following this 
line of thought one could further use the 8 C (M, z) predicted by the 
extended spherical collapse model as the correct moving barrier. 



> 




M h [M Sjn /h] 

Figure 4. Ratio between the PS multiplicity function evaluated with 
S c extended and 5 C acdm (upper panel) and ratio between the PS mul- 
tiplicity function with <5 C extended and the ST multiplicity function with 
S c , ACDM - Line styles and colours are as in 



Therefore by dire ct comparison of expressions (2), (4) and (5) in 



Sheth et 



c by dtrec 
all feOOll) 



we substitute S c , extended directly into their ex- 
pression (2) (equivalent to the PS multiplicity function). Results are 
shown in Fig.f4]for the three redshifts of Fig. [3] 

In the upper panel we show the ratio of the PS multiplicity 
function evaluated with the S c , extended and 5 c ,acdm linear over- 
density parameters while in the lower panel we compute the ratio 
of the PS multiplicity function evaluated with the Sc, extended linear 
overdensity parameter with the ST multiplicity function evaluated 
with the <5 c ,acdm linear overdensity parameter. As it is evident, 
when we compute the ratio using the same multiplicity function, we 
obtain a very similar behaviour as in Fig. [3] even if quantitatively 
slightly different. When we instead identify rotation with ellipticity 
and use directly the PS multiplicity function and compare it with 
the ST multiplicity function, we observe a totally different behav- 
ior. This can be explained not only in terms of different overdensity 
parameters, but also remembering that the PS mass function pre- 
dicts more (less) objects at low (high) masses with respect to the 
ST parametrization. Our results also show that if we want to repro- 
duce the results of the ST mass function using the extended spheri- 
cal collapse model, we need to modify the moving barrier and find 
a new parametrization for the multiplicity function. This goes be- 
yond the purpose of this work, therefore, taking into account these 
caveats we will assume the ST mass function as the correct one. 
In Fig. [5] we show the ratios, at the three different redshifts 
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considered, between the different dark energy models here consid- 
ered and the extended ACDM model, where the term a 2 — u> 2 is 
included. 

As we can notice, the (generalised) Chaplygin gas shows a 
huge suppression of structures at all redshifts, making therefore 
this model ruled out in the extended SCM. All quintessence models 
have a lack of high mass objects. While this is not severe at all for 
the INV2, 2EXP and CNR model, all the others have a suppression 
of several orders of magnitude, increasing with the increase on the 
redshift. The most suppressed model is the AS model, that shows 
also the most different S c from the extended ACDM model. 
Regarding the phantom models, differences are at most of a fac- 
tor of 4-5. While the model with w = —2/3 shows a decrease 
in structures, the phantom models show an increase. Differences 
are significant in general only for high masses M > 10 14 Mq/Ii 
while for the model with w — —2/3, they are evident already at 
M « 10 13 for z = 1. 

We also notice that a small amount of curvature has a very little 
effect on the number of objects, as differences are of the order of 
few percent even for cluster scales. 

Out results can be easily interpreted in terms of the differ- 
ent matter power spectrum normalizations. The Chaplygin gas has 
an extremely low normalization (erg = 0.066) making therefore 
very unlucky that structures could form in such a universe. Phan- 
tom models instead show an higher normalization, making there- 
fore easier to have high mass objects. Moreover one has to take 
into account that now the linear overdensity parameter 8 C is modi- 
fied and very strong differences will reflect in the differential mass 
function (see for example Fig.|2j- 

A direct co mparison between the results in this work and the 
ones described in IPace et al.l j2O10h cannot be made for several rea- 
sons. In particular, the power spectra here were normalised in or- 
der to h ave the same amp litude of fluctuations at the CMB epoch 
while in IPace et al.l J201C|) a different normalisation was adopted. 
There the power spectrum was normalised in order to have nearly 
the same mass function at z = 0. This implies that effects of dark 
energy will be important only at high redshifts while in our case we 
see substantial differences already at low redshifts as expected. In 
addition, here we just limit ourselves to the study of the differential 
mass function and we do not investigate the cumulative number of 
halos. This is because we do not want to have our results affected 
b y volume effects. N ote also that due to the choice of normalization 
in IPace etal.ll2O10l) . the quintessence models will predict more ob- 
jects than the ACDM at high redshifts (not shown here). 



4 CONCLUSIONS 

In this work, we study the impact of the term a 2 — uj 2 on the spher- 
ical collapse parameters, namely the linear overdensity parameter 
S c and the virial overdensity parameter Ay and how this reflects on 
the number of objects via the mass function formalism for a broad 
cl ass of dark energ y models, already studied in the spherical limit, 
bv lPace et al]d2010l) . 

We assume that only the dark matter component is cluster- 
ing and that dark energy is only at the background level, therefore 
affecting only the time evolution of the Universe. Doing this, we 
implicitly assumed that eventual perturbations in the dark energy 
component can be neglected. 

We showed that the non-linear term considered opposes to the 
collapse and this is reflected by higher values of the linear overden- 
sity parameter with respect to the spherical case. Modifications are 



quite substantial, of the order also of 40% for the ACDM model. In 
general the effect of dark energy is to lower the value of <5 C with re- 
spect to the ACDM model and we see that this is also the case also 
in the extended SCM. Despite the values of the linear overdensity 
parameter are higher now than in the non-rotating case, we see that 
in general dark energy still lowers its value. This is the case only if 
w > — 1, while for phantom models the super- negative equation of 
state slows down the collapse. 

In order to appreciate the interplay between the different dark 
energy equation-of-state parameter and the power spectrum norm- 
lization, it is interesting to compare models having approximately 
the same power spectrum normalization. For example, this is the 
case for the models INV1, CPL and Casimir (see Table [T] for the 
exact values). Having very similar power spectra normalizations, 
from Fig. [5] we see that the ratio with the ACDM mass function 
gives in general very similar results. Models INV1 and CPL are 
very similar to each other and in general similar to the Casimir 
model. Differences between these models can be seen more clearly 
in the evolution of S c in Fig. [2] This shows that often differences 
in the models are hidden by the power spectrum normalization. In 
general, taking into account this caveat, for other models differ- 
ences in the differential mass function are due to the influence of the 
dark energy component (see also comments at the end of Sect. [3}. 

As expected, such differences reflect in the number of ob- 
jects. Since with respect to the standard case the only quantity to 
be changed is <5 C , we can easily study the impact of a rotation term 
on structure formation. We show this in Fig. [5] The term a 2 — ui 2 
suppresses, as expected, the high mass tail of the mass function, 
since rare events are more sensitive to the background cosmology 
and to the collapse process. In general low masses objects are not 
severely affected by rotation, but a noteworthy counterexample is 
given by the (generalized) Chaplygin gas and AS model where we 
observe a suppression in the number of objects already of several 
orders of magnitude for galactic masses. 

We conclude therefore that the term a' 2 — lu' 2 has a strong 
impact on structure formation and that it is worth to investigate 
different parametrizations for the additional term. 
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